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Xf^ , Abstract 

ff^ ■ In this paper we present a two-field inflation model, which distin- 

Tjlj- . guishes itself with a non-canonical kinetic lagrangian and comes from 

CD \ the large volume approach to the moduli stabilization in flux com- 

pactification of type IIB superstring on a Calabi-Yau orientifold of 
/),(i.2) ^ /i(i.i) > 4. The Kahler moduli are classifled as volume modu- 
lus, heavy moduli and two light moduli. The axion-dilaton, complex 
/\ ' structure moduli and all heavy Kahler moduli including the volume 

C^ ■ modulus are frozen by nonperturbatively corrected flux superpotential 

and the a'-corrected Kahler potential in the large volume limit. The 
minimum of the scalar potential at which the heavy moduli are sta- 
bilized provides the dominant potential energy for the survived light 
Kahler moduli. We consider a simplified case where the axionic com- 
ponents in the light Kahler moduli are further stabilized at the poten- 
tial minimum and only the geometrical components are taken as the 
scalar fields to drive an assisted-like infiation. For a certain range of 
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moduli stabilization parameters and inflation initial conditions, we ob- 
tain a nearly flat power spectrum of the curvature perturbation, with 
Ug ~ 0.96 at Hubble-exit, and an inflationary energy scale of 3 x 10^^ 
GeV. In our model, significant correlation exists between the curva- 
ture and isocurvature perturbations on super-Hubble scales so that 
at the end of inflation a great deal of the curvature power spectrum 
originates from this correlation. 

1 Introduction 

The development of top-down cosmology has encountered a research hot on 
studying the superstring landscape for realizing the cosmological inflation. 
This involves the search of metastable de Sitter-like string vacua at which the 
ten-dimensional spacetime behaves effectively as having only four dimensions, 
with six spatial extra dimensions compactified on a Calabi-Yau threefold and 
almost all of the compactification moduli including the Calabi-Yau volume 
stabilized. The advent of KKLT mechanism[T] for moduli stabilization has 
opened lots of possibilities to engineer inflation in string theory, using either 
the hght closed string modulijH El IH [S] or open string modulij^ related 
to the position of a mobile D-brane in wrapped geometry as the relevant 
scalar flelds. Different stringy realization of cosmological inflation may not be 
mutually incompatible. They may arise in different regions of the landscape, 
among which we anticipate that there is at least one possibility that realizes 
the inflationary evolution of our Universe. The different inflation scenarios 
may also turn into one another, increasing the overall probability of inflation 
from a single mechanism [5]. 

In a reasonable inflationary scenario, the inflaton fleld has to be ensured 
to roll its potential very slowly so that there are enough efoldings for ac- 
celeration. It is necessary that, during the inflation epoch, the potential 
energy curve of the involved scalars must be sufficiently flat and their effec- 
tive masses have to be much smaller than the Hubble parameter, m'^ -C H^. 
Scalar fields such as the brane inter-distances in a general D-brane inflation 
model are conformally coupled to gravity, which enables them to acquire 
squared masses m^ ~ j^R ~ H^ (//-problem) . To overcome this r^-problem 
and have a slow-roll inflation in such a scenario, severe and sometimes un- 
natural flne-tunings are indispensable [6]. Alternatively, no severe ?7-problem 
exists in the closed string moduli inflation scenario[2],[3llll[5] where the infla- 



tion is thought of to be driven by some hght closed string moduh associated 
with 4-cycle volumes of the compact Calabi-Yau threefolds, particularly in 
the large volume limit [7] |8]. 

The first proposed closed string moduli inflation scenarios in Type IIB 
superstring theory are the racetrack inflation model [2] and the better race- 
track inflation model[3], where the heavy moduli are stabilized at a super- 
symmetric anti-de Sitter F-term potential minimum and the potential is up- 
lifted to have a de Sitter vacuum for inflation in an uncontrollable manner 
unless some additional ingredients are added to break the supersymmetry 
spontaneously [9l [TOl [HI [12]. In the subsequent Kahler inflation model[l] 
and Roulette inflation model[5|, the scalar potential for the moduli flelds is 
calculated in the large volume limit, where the leading a'-correction to the 
Kahler potential has been taken into account so that the heavy moduli are 
stabilized at a supersymmetric broken anti-de Sitter F-term potential mini- 
mum. Because this anti-de Sitter minimum is apart from supersymmetry, it 
can be uplifted to a de Sitter vacuum by taking into account some D-term 
corrections to the scalar potential[T31 [T3]. Moreover, after heavy moduli are 
stabilized in the large volume limit, the scalar potential for the survived light 
Kahler moduli are exponentially flat, making such approaches very promising 
for engineering assisted-like inflation [15]. 

In string theory cosmology scenario, the early evolution of our universe is 
essentially drived by multiple scalar flelds. When superstring is compactifled 
on a Calabi-Yau threefold to yield a four-dimensional effective M = 1 super- 
gravity theory, lots of moduli flelds emerge. It has been conjectured, however, 
that increasing the number of Kahler moduli flelds does probably make the 
inflation be easier to achieve [3]. Moreover, the presence of multiple scalar 
flelds during inflation can lead to quite different inflationary dynamics such 
as the detectable non-Gaussianity in primordial density perturbations on 
super-Hubble scales and residual isocurvature fluctuations after inflation [T6] 
that might appear unnatural in a single fleld model. 

Our aim in the present paper is to build a two-fleld inflation model from 
the large volume limit of type IIB superstring with 3-form fluxes compact- 
ifled on a Calabi-Yau orientifold of /i*^^'^) > h^^'^^ > 4 and discuss its impli- 
cations on COBE spectrum of the linear curvature perturbation. For that 
we flrst investigate the possibility of realizing better racetrack inflation in 
the large volume approach. We divide the Kahler moduli into volume mod- 
ulus, h^^'^^ — 3 heavy moduli and two light moduli. The scalar potential is 
dominated by its F-term contribution, and is corrected by a volume-modulus 



dependent uplifting. The axion-dilaton, complex structure moduli and all 
heavy Kahler moduli including the volume modulus are stabilized, in the 
large volume limit, to a de Sitter-like minimum of the potential, whereas 
the two light Kahler moduli remain dynamical during the moduli stabiliza- 
tion process. These two light Kahler moduli are supposed to be the scalar 
fields driving the inflationary expansion of the 4-dimensional universe. The 
de Sitter minimum of the dominant potential does not only stabilize the 
heavy moduli by giving them masses, it governs also the evolution of two 
survived light moduli. The exponential dependence of the nonperturbatively 
corrected superpotential upon the Kahler moduli guarantees the flatness of 
the scalar potential for two light Kahler moduli. If the axionic components 
of the light moduli are dynamical, this is a large volume version of the better 
racetrack model f3] in which four scalar fields are involved. In the present pa- 
per we choose to stabilize the axionic components at the potential minimum, 
a two-field assisted-like model[15] results in instead. In some sense, this is an 
extension of the single field inflationary model|4| to two-field case, to which 
we calculate the linear perturbations in a strict multi-field approach. By 
assigning some suitable model dependent parameters, we obtain the COBE 
normalization favored power spectrum of the linear curvature perturbation, 
and in particular, Ug ~ 0.96 at the horizon exit. Dependent upon the tra- 
jectory in phase space and thus the initial conditions we choose, the isocur- 
vature perturbation in this model could decay very slowly for most of efolds 
on superhorizon scales. We have also observed that the power spectrum 
of curvature perturbation on the super-Hubble scales has some remarkable 
departure from that of adiabatic perturbation at the end of inflation, im- 
plying that the curvature fluctuation after inflation does inevitably include 
significant contributions from isocurvature perturbation. 

The paper is organized as follows. In the second section, the theoretical 
setup is introduced. We review briefly the large volume approach to the 
moduli stabilization developed by Quevedo and his collaborators in [H [8] 
and establish our own two-field model. The third section is devoted to the 
investigation of possible inflationary dynamics in our model. We first present 
the equations of motion for the background fields and the linear perturba- 
tions, and then calculate the power spectra of these perturbations and their 
correlation by numerical integration. The obtained curvature power spec- 
trum is nearly flat during most inflation and is characterized by a spectral 
index Ug ^ 0.96 at the claimed Hubble crossing. The flnal section is our 
conclusions. 



2 Inflationary landscape in large volume ap- 
proach 

Inflation in string theory is always associated with the resolvability of the 
moduli stabilization because the potential used to stabilize the heavy moduli 
might yield unacceptably large masses for the inflation driversjO]. In type 
IIB superstring compactified on a Calabi-Yau orientifold, the closed string 
moduli consist of the complex structure moduli, dilaton field and the Kahler 
moduli. All complex structure moduli including the dilaton field could be- 
long to the heavy moduli and be stabilized at a minimum of the effective po- 
tential if we introduce an imaginary self-dual 3-form flux into the orientifold 
construction [TTIIIH]. The stabilization of the Kahler moduli requires, accord- 
ing to KKLT mechanism[l], introducing the nonperturbative contributions 
of the Euclid D3-branes or some wrapped D7-branes to the superpotential. 
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where Tj = Tj + iOi with Tj the i-th 4-cycle volume and 6i the corresponding 
axions. The coefficients Ai represent threshold corrections and are indepen- 
dent of the Kahler moduli. 

The Kahler potentials arising from type IIB superstring are of no-scale 
or approximately no-scale [ 
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K = K,s-2\n{V+^] (2) 
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K'^diKd-jK - 3 = ^^, ^ ^ (3) 



in the large volume limit V — > oo, where f = — ^^ '^^ i,L , C(3) ~ 1.2 and 

x{M) = 2{h^^'^^ — h^^'"^^). In ([2]) the a'-corrections have been included. Com- 
bining with the superpotential given in Eq.([T]), we see that the Kahler moduli 
do only appear exponentially in the F-term scalar potential, 

Vf = e^'iK'W.KDjK-SlWl') 



^Throughout we take the units Mp = I/vSttG = 1. 



where DiK = diKW + diW . The uphfting correction to the scalar potential, 
which has several possible sources and is necessary to fine-tune the cosmo- 
logical constant, depends on moduli only through the overall volume, 

Vupiift ~ Ye (^) 

where | < f? < 3. For concreteness, we assume ^ = 2 in the present paper. 
It is remarkable that, in the presence of several Kahler moduli, the variation 
of scalar potential along each Ti direction is in general uncorrelated with 
its magnitude. Besides, the potential is exponentially fiat along some Tj 
directions as long as these Kahler moduli are sufficient large. Therefore, the 
large Kahler moduli have much chance to become the scalar fields to drive 
the early infiationary evolution of our Universe. 

The Calabi-Yau orientifolds of Swiss-cheese type have been shown to be 
very useful in string infiation engineering, whose volume can be formulated 
into a simplified form as follows [H [5], 
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Here the volume of the i-th 4-cycle is denoted by Xj = SfJTj, among which 
Ti controls the overall volume and T2, ■ ■ ■ , r„ [n = h'^^'^^] are blow-ups whose 
only non-vanishing triple intersections are with themselves. We stabilize the 
dilaton and complex structure moduli with fiuxes, following the plausible 
KKLT procedure[lj. After that, the superpotential ([T]) is reduced to be only 
Kahler moduli dependent, W = Wq + Yl^=i Aie~"''-^\ where Cj = 2'k/N and 
for simplicity Wq is assumed to be real. 

We establish our formalism in the large volume limit, where V — > cxd, 
Ti ^ V^/^ -^ 00, ajTj ~ InV for z = 2, 3, ■ ■ ■ , ra but Ai ^ Aj for i = 
2, 3, ■ ■ ■ , k and j = k + 1, k + 2,- ■ ■ , n. Because A^ in superpotential ([T]) 
are directly proportional to the squared masses of the corresponding Kahler 
moduli, Ti are referred to be the heavy moduli when i = 2, 3, ■ ■ ■ , k and light 
moduli when i = k + 1, k + 2,- ■ ■ , n, respectively. In this limit, e^ ~ ^f^"^, 
the superpotential is approximately independent of Ti which we will call the 
volume modulus, 

n 

W^Wo + J2 Ae""'^^ (7) 

1=2 



Up to the magnitude of order 1/V^, we can approximately write the scalar 
potential as l^ ~ Vdom + Vcorr, where the dominant part Vdom consists of 
the contributions of the volume modulus and the heavy moduli, and of the 
necessary uplifting potential, 
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The contributions of the light Kahler moduli to the scalar potential do only 
appear in correction term Vcorr 1 

y^'- t '-^^M^^-'"-+ ± '-^^^-cosiaA, (9) 

i=k+l * ° i=k+l 

In Eq.®, we have parameterized the uplifting potential as VupUft = ^yr-, 
with 7 a positive parameter. There are terms not included in Eqs.® and 
(111), which are sub leading. 

In the large volume approach [3, [8], the moduli stabilization problem of 
the Kahler moduli is solved in the following two-step procedure [20]. First, we 
stabilize the axions 6i (z = 2, 3, ■ ■ ■ , fc) to the potential minimum by setting 
ttj^j = 71. Second, we can find an approximate minimum of the potential 
by letting Vdom be flat along the directions of V and the heavy moduli Tj 
{i = 2,3, ■ ■ ■ ,k) in moduli space. The minimum of the scalar potential turns 
out approximately to be. 
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where V/ that stands for the stabilized Calabi-Yau volume is the solution of 
the algebraic equation 

^^_9Sy_^f^W^EIdl2,.V^(2c,.in (11) 
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with 

Ci = In 






The heavy moduh will be stabilized at aiTimin ~ InV — Cj. Phenomenolog- 
ically, we expect that Vmin provides the dominant potential energy for the 
inflationary evolution of the survived light Kahler moduli. The COBE nor- 
malization of density perturbations at the Hubble exit 6h ~ 1-92 x 10~^ then 
demands(Kii„/ei)^ ~ 0.004 ~ 5.4 x lO^*^ GeV at the Hubble exit, where ei is 
the first slow-roll parameter. The range of ei at the horizon exit is ei ~ 10~^^ 
for typical values of the microscopic parameters. Thus the inflationary energy 
scale should be rather low in such models, Vmin ~ 10~^^[4]. This requirement 
can be easily satisfied by making some fine-tuning on the involved stringy pa- 
rameters and, in fact, there exists a landscape in parameter space in favor of 
phenomenology and cosmology. One possibility that we will focus on in the 
next section is given by an orientifold model of h^^'^^ = 4, h^^'"^^ ~ 100, with 
geometrical parameters ao = 1/9a/2 and A2 = A3 = A4 = l[5j. The complex 
structure moduli including dilaton are supposed to be fixed by imaginary 
self-dual 3-form fluxes in such a manner that after these moduli are stabi- 
lized we have Wo = 100, gs ~ 0.132 and consequently C, ~ 10. The other 
microscopic parameters in the model are designated as, 

A2 = 1, A3 = A4 ^ 1.075 X 10~^ aa = 03 = 04 = 27r/300, 
7 ^ 1.026 X 10-^ 

For such a model, numerical calculation tells us that V/ ~ 2.8 x lOij and 
Vmin ~ 4.035 X lO^^Tj. The volume modulus ri and heavy modulus T2 are 
found to be frozen at nj ^ 4.28 x 10^ and 02X2,/ ^ 13.81 ^ 0.71 In V/ 
respectively, as expected. 

The large volume approach to heavy Kahler moduli stabilization naturally 
provides a platform for having a stringy inflation scenario, where the survived 
light Kahler moduli, say Tj for i = k + 1, ■ ■ ■ , n {k > 2) play the roles of the 
multifield inflation drivers. The effective potential for these light moduli is 

V « V„.,.+ ± M^e--+ ± l^^lj^e-"" cos(a.».) (12) 

and these scalar fields are assumed to be initially far from the minimum of 
potential V (Please distinguish Vmin from the minimum of V. The latter is 



^That is, Vf « 2.8 x 10^ /f , where h = 2TrVa'. 

•^Correspondingly, the density of the vacuum energy during inflation is of the scale 

M « VJ^j„ w 3.075 X 10^** GeV. This setup is consistent with the upper bound V^ < 
3.8 X lO^s GeV from the present WMAP dataset[2T]. 



expected to be zero) . Because in Eq. (IT^ Vmm is one order of the fixed Calabi- 
Yau volume in the magnitude greater than the other scalar field dependent 
terms, the fixed volume V/ and all of the stabilized heavy Kahler moduli 
as well as the associated axions will remain frozen during inflatioio. These 
frozen heavy moduli will not create the harmful contributions to the masses 
of the light scalar fields. 

The form of potential flT^ implies that there is a stringy landscape to 
realize the cosmological inflation in IIB flux compactification at the large vol- 
ume limit. If we consider the type IIB compactification on the Calabi-Yau 
orientifolds with h^^'^^ = n = 3, for example, we will get Kahler modu- 
lus inflation |3] when the axion 63 is stabilized at the potential minimum or 
roulette inflation [S] when 63 is dynamical. In this paper, we consider the IIB 
compactification on a Calabi-Yau orientifold of h^^'^' = n > A, with almost 
all Kahler moduli (including V, ri and Ti = Ti + i9i ioi i = 2,3, ■ ■ ■ ,n — 2) 
have been stabilized. The only exception is for the two light moduli T„_i = 
Tn-i + ^^n-1 and Tn = Tn + iOn which correspond to four scalar fields. The 
dynamics of this four-field system is governed by an effective potential 
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T/ ~ T/ . _i_ V 2iffi£li2_yZip- 

n '=''-' (13) 

+ E '"'y"' e-"'"'cos(a,g,) 

i=n—l f 

Recall that the effective superpotential for these fields is a sum of two ex- 
ponential terms, the potential flT^ might be a large volume version of the 
better racetrack inflation model[3j. 

The detailed investigation to this refined better racetrack model is sched- 
uled to be published in the near future[22]. In the present paper, we do 
what Conlon and Quevedo have done by fixing the axions to the potential 
minimum and ignoring the double exponentials, but regard both light Kahler 
moduli T^ and T4 as the relevant scalar fields (Here we take n = 4 for sim- 
plicity). This will yield an assisted-like extension of the single field Kahler 
modulus inflation model. The potential energy for such a two-field system 



""^The unique exception is the axion associated with the volume modulus, 9i. It is free 
of stabilization in the large volume approach. However, this is harmless. 



The kinetic Lagrangian is given by £^k = — X]j j=sKijd^Tid'^Tj. Explicitly, 
Let 



/4aoAi+2 3/4 /.-x 

^i = \l gy ^i+2 (15) 



with « = 1, 2, we get, 
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^K = --d'^ipidf.ipi - -d^'ip2d^(p2 - -ViV2d^ipid^,Lp2 (16) 

The potential (1141) can be recast as 

V{ipi, ip2) = Vrran + Vi'ilj{ipi)e-''''^^''^^ + V2ij{ip2)e~^''^^'''''^ (17) 

in terms of the newly defined fields. In Eg. (1 171) we have defined 4'{'^i) = </?/, 
Pi = 0^+2( 4^^/ ) 3 and Vi = _ 4ft '+2^0 fQj, convenience. 



3 Two-field infiation model with non-canonical 
lagrangian 

In this section, we intend to study the probable implications of the above 
two-field model on cosmology. The model emerges from the large volume 
approach to the moduli stabilization of Type IIB superstring fiux compacti- 
fication and is described by the effective action 

s= J d^xV^\^-n - ^d^'^id^^i - \d^^2d,^2 

9 ] 

- ■^^iV>2d^(pidf,(p2 - V{ipi, 1P2) (18) 

o 

with V{(pi,ip2) given in Eq. flTTl) . Similar to the assisted infiation model[T5]. 
the potential V{ipi, ip2) consists of the sum of two exponentials of the scalar 
fields. The model can be thought of as a simplified version of the large volume 
better racetrack model in which two axionic fields have been stabilized at 
the potential minimum. In action fITS]) . G stands for the determinant of 

10 



4- dimensional metric tensor 0^,^, whose Ricci scalar curvature is denoted 
by TZ. From the assigned stringy parameters for the model, the involved 
phenomenological parameters take values of Vmm ~ 4.035 x 10~^^, Vi = 
V2 ^ -2.212 X 10-1^ and /3i = /^a ^ 40337.578 respectively. The potential 
has a Minkowski minimum, V{(fi,ip2) ~ 0, at yji = 992 ~ 3.51 x 10~ij, 
however, it is very fiat far away from this minimum. 



V(^, X) 




0.000 



Figure 1: The potential surface of our two-field model defined in Eqs. ()20p and 
(pT|) where the values V^in ~ 4.035 x IQ-^^ 14 = ^2 ~ -2.212 x IQ-^^ and 
/?! = /32 ~ 40337.578 have been assigned to the parameters. The potential has 
a Minkowski minimum at (/? ~ 4.97 x 10~^ and x = ^^^ is very flat along the 
99-direction far away from this minimum. 



The crossing term in the kinetic lagrangian brings on much unnecessary 
inconvenience. To overcome it, we redefine the independent scalar fields by 

^1 = {^ + X)/V2, ^2 = {^-x)/V2. (19) 

The action becomes 

-7^ - cos^ a d^ifd^if - svo? a d^x^^tX - '^(v', x) 



S 



d^xV^ 



^Correspondingly, the light Kahler moduli at this Minkowski minimum take values 
Ta = r4 w 47.75. 
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(20) 
where 

and cos2q; = |(y9^ — x^). Because in large volume limit the magnitudes 
of both Lf and x could be much smaller than one, the introduction of the 
field-dependent auxiliary quantity «(</?, x) by its cosine value is reasonable. 
The appearance of a{ip,x) iii the kinetic lagrangian implies that we get a 
nontrivial diagonal metric in field space, which does not coincide with the 
known non-canonical two-field kinetic Lagrangian studied in [23], and is ex- 
pected to bring some new features in cosmological application. In terms of 
newly defined scalars, the Minkowski minimum of the potential (12T!) occurs 
at c/9 ~ 4.97 X 10~^ and x = 0. Departure from this minimum is a plateau 
mainly along the yj-direction which might be adequate for inflation. 

We begin with the equations of motion of the homogeneous background 
fields. Consider a spatially flat Robertson- Walker spacetime 

dsl = -dt^ + a{tfdx^ (22) 

Here t is the cosmic time. It is much convenient to formulate the equations 
of motion of the scalar fields in terms of the efold time n = In (a(t)/aj„j), 
by which the equations of motion of the scalar fields are decoupled from 
the metric evolution. For the sake of convenience, we define a homogeneous 
infiaton field a(t) through its velocity. 



a = ^2cos^a<^^ + 2sm^ax^ (23) 

In Eq. (l23p and hereinafter, a dot stands for a derivative with respect to the 
efold time n. The equations of motion of the scale factor in metric (!22l) and 
the homogeneous background scalar fields are found to be, 

I - -¥ p«' 

(^ — ta:^ ad^a{ifP' + X^) — 2 tail ad^aipx . sec^ ad^pV 

3 - a72 +^= W ^ ^ 

X + cot ad^a^ifP' + x^) + 2 cot a^^paipx ■ csc^ c^^x^ 

3 - ^2/2 "^ ^ " 2V ^ ^ 

12 



where V is the abbreviation of the potential V{(p,x), H = da/adt is the 
Hubble parameter but H = dH/dn. Inflation occurs for ^ > 0. Hence, 
having an inflation driven by the above two scalar flelds requires 

ei = -^ = crV2 < 1 (28) 

ei is the so-called flrst slow-roll parameteio. Is there an inflationary epoch 
with ei < 1 for our two-fleld model? To have a deflnite answer to this 
question, we have to numerically solve Eqs. (12411271) under appropriate initial 
conditions. In fact, it is sufficient to integrate Eqs. (1261) and (l27j) only, since 
these scalar field equations are decoupled from the metric evolution. 

From the Cauchy theorem, the solution to Eqs. (1261) and (I27p is unique 
provided the initial fields and initial field velocities are given at some initial 
instant n = Uini- However, just as pointed out by Ringeval[2^, the attractor 
behavior induced by the friction terms erase any effect associated with the 
initial field velocities after a few efolds. The integration of Eqs. (1261) and (1271) 
depends essentially upon the initial field values only. The attractor behavior 
does also ensure the stability of forward numerical integration schemes, so 
we use a Runge-Kutta integration method of order four. Numerically, we 
take Uini ^ -4.605ld3 and (pini ^ 5.67 x lO'^. It follows from Eq.([T9]) that, 
if Xini vanished, fi^ini = ^2,ini, the evolutions of ipi and (p2 would exactly 
be the same, our multifield model would effectively reduce into a single field 
model. To have a small but significant deviation from such an indifferent 
situation, we take Xmi ~ —1.13 x 10~^. Such a choice for initial fields ensures 
them far away from their possible values at the Minkowski minimum, and 
correspondingly the light Kahler moduli are initially set to be t^ ~ 1224.14 
and T4 ^ 1230.67, respectively. The initial field velocities are chosen on the 
attractor by setting 

^ini = -d^ In l^|^,„„ x,„, X^ni = "^x !» ^Iv.n., X^n^■ (29) 

Figure 121 gives our numerical solution to the classical trajectory of the scalar 
fields for -4.60516 <n< 56.83847. During most of the efolds, -4.60516 < 
n < 55.85, the scalars roll very slowly and the trajectory is almost in the 
{(p — x)-direction. Being a two-field model means that the trajectory is 



^The slow-roll approximation is in general figured in terms of numerous Hubble flow 
functions [24j. among which the first two are ci = — § and 62 — j^- 
^See Eq. ([Tf|) below for explanation. 
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not a straight line and in fact a sharp turn on trajectory occurs roughly 
at n K. 55.85. Within the last efold time interval the trajectory becomes 
another straight line along the (</? + x)-direction. Plugging this solution for 
(p{n) and x(^) i^to Eqs. flMj) and fl25l) determines the evolution of Hubble 
parameter and the first slow-roll parameter ei. Our choice of the initial 
conditions ensures both ei and €2 much smaller than one for most of the 
efolds. In particular, ei ^ 4.82 x 10~^^ and 62 ~ 0.036 at the Hubble exit 
(n ~ 2). There is indeed an inflationary epoch in our model during which 
the Hubble parameter hardly changes. The potential V is very fiat along the 
classical trajectory and, just as in the case of single field case[l], it is almost 
a constant, V ~ 4.035 x 10^^^, for most of inflation. Only within the final 
fractions of the last efold does a steep decay of potential take place. When 
n impends over 56.83847, the potential sharply drops off from the preceding 
constant to a metastable minimum with V ~ 2.725 x 10~^^, where ei ~ 0.97 
and the inflation is close to its end. 

The apodictic deviation of trajectory from an entire straight line shows 
that what we consider is a genuine multifield model. Provided that the 
isocurvature perturbation decays fast enough on the super-Hubble scales, it 
will not influence the super-Hubble curvature perturbation heavily even if 
the trajectory in field space deviates deadly from the straight-line one of a 
single field model. The steep decay of the isocurvature perturbations takes 
place in many supergravity or string theory inspired multifield models, and 
in particular in the string theory inspired roulette model [5]. The infiation 
in these models is effectively driven by a single scalar field where one can 
use the single field approximation to account for the super-horizon power 
spectrum of the curvature perturbation. However, the same is not true for 
the present situation. In our two-field model, the decay of the isocurvature 
perturbation on super-Hubble scales is not fast and the spectrum of the 
curvature perturbation in the vicinity of the Hubble crossing is remarkably 
deviated by that at the end of infiation, where the correlation between the 
curvature and isocurvature modes becomes sufficiently strong so that the final 
curvature perturbation is heavily sourced by the isocurvature perturbation. 

We are going to confirm the above conclusion for the two-field model 
in consideration through numerical calculations of the power spectra of the 
linear perturbations and their correlation. Due to the fact that at the Hubble- 
exit ei ^ 4.82 x 10~^^, there is no important tensor fiuctuations in our model. 
This might be a common feature of the stringy infiation models where the 
inflation is driven by some closed string moduli[21 El IH [S]- So we consider 
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only the scalar perturbations. Because the matter is composed of two scalar 
fields, the stress tensor is diagonal and the perturbed metric is of the form 



dsi 



[1 + 2^(t, x)] dt^ + a{tf [1 - 2^(t, x)] rfx^ 



(30) 



in the longitudinal gauge. The scalar fields are decomposed into their homo- 
geneous backgrounds plus the anisotropic perturbations: 



ip{t, x) = if{t) + 5ip{t, x), x(^, x) = x{t) + ^X(^, x). 



(31) 



To facilitate solving the evolution equations of these cosmological perturba- 
tions, we can alternatively decompose them into the instantaneous curvature 
and isocurvature components, relying on the fact that the isocurvature modes 
do only source the curvature perturbations approximately |25l l26j . The cur- 
vature and isocurvature perturbations are parallel with and orthogonal to 
the trajectory of the homogeneous inflaton field cr(t), respectively. To define 
them we first introduce a time dependent angle 6', 



cos6' 



v2cosa — , sin0 = v2sina;— . 



a a 

with which we can recast the background equations (l26l) and (1271) as 

a + (3-aV2)(T+^ = 



e + 



V, 



a 



H^a V2 



m 



l^cos 6 sec aa-^ + sin 6 esc aa^p) = 



(32) 

(33) 
(34) 



The curvature perturbation 6a{t, x) and isocurvature perturbation 6s{t, x) 
are defined through a linear transformation in field space 



5cr(t,x) 

6s(t,x.) 



V2 



cos a cos 6 sin a sin 6 
— cos a sin 6 sin a cos 6 



5v3(t,x) 

5x(^,x) 



(35) 



From Eq. (l5B]) we have 6s = sm2a{ip6x — X^^^jb. The isocurvature pertur- 
bations turn out to be the entropy perturbations. 

It is worthwhile to stress that, among the linear scalar perturbations 
\E'(t, x), (5cr(t, x) and (5s(t, x), only two of them are independent. The in- 
dependent scalar perturbations can also be chosen as the gauge invariant 
Mukhanov-Sasaki variable 



Q^(t, x) = b(y{t, x) + (T^(t, x) 



(36) 



15 



and (5s(t,x). In terms of these gauge invariant perturbations, the perturbed 
Klein-Gordon equations and perturbed Einstein equations [21] can be unified 
into two coupled second order differential equations: 

g. + (3 - dV2)g. + j^6s - ^^v^g. + ^g. + ^6s = m 



■ 2V ■ 
5s+i3-&y2)5s-y^Q, 



H^ 



a 



_L_V^,.+ ||Q, + |l,. = 0(38) 



where 



V^^ + H'\3-ay2)a'- — {^ 



2&V„ 



+ 



V2K 
sin^ 2a 

V2Vs 
sin^ 2a 



C 



^2^K 2K^ 



a 
2\/2K 



(cos 29 — cos 2q;) (cos 9 sin aa^ + sin 9 cos aa^^,) 

cos 9 cos a(cos 26' + cos 2a — 2)a^ 

— sin 6* sin a (cos 26^ + cos 2a + 2) a^ (39) 

+ 2Ks 



+ 



C. 

c. 



sin^ 2a 
2V2K 
sin^ 2a 



(cos 26* — cos 2a) (cos 9 cos aa^, — sin ^ sin aa^) 
(cos 26^ + cos 2a) (cos 9 sin aa;^ + sin 9 cos aa^|40) 



a H^a' 



-6— - 



+ Kcr 



(41) 



V2Vs 



sm 2a 



+ 



sin^ 2a 

V2K 

sin^ 2a 



cos 26* + cos 2a) (cos 9 sin aa^ — sin 9 cos aa;^) 

cos 9 sin a(2 — cos 29 + cos 2a)a<^ 

— sin^cosa(2 + cos26' — cos2a)a;^ (42) 



As usually done, we will work with the Fourier components of the perturba- 
tions, Q\i{n) and 5sk(^), with k a given comoving wave number. For sake 
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of convenience, in Eqs. fl37ti^2l) we have defined various derivatives of the po- 
tential with respect to the curvature and isocurvature directions, which are 
associated with the derivatives of the potential with respect to the original 
scalars if and x through the linear transformations 



Vs 



r 



and 



K. Vr. 






T 



V^ 
K 



V.o.n K 



'ipip 



vx 



K.V K 



V'X 



XX 



r^ 



The transformation matrix T reads, 



^^71 



cos 9 sec a sin 6 esc a 
- sin 9 sec a cos ^ esc a 



(43) 



(44) 



(45) 



To study the generation of perturbations from the vacuum fluctuations, 
we have to solve the closed set of Eqs. (157|) and (1551) under appropriate initial 
conditions. Though the curvature and isocurvature fluctuations are finally 
coupled to each other beyond the Hubble crossing, they might be statisti- 
cally independent deeply inside the Hubble radius. The initial efold time 
Tiini at which the perturbations are thought of as mutually independent must 
be much less than the Hubble crossing instant n^, for Fourier mode pertur- 
bations of wave number k, k < a{n^.)H{n^). The value of n^, depends on 
the efold number during which the universe reheated before the starting of 
the radiation era[24j and it is known [2 7j that typically 40 < A'^* < 60 for 
N^ = Uend — n*- Without knowledge of details of the reheating mechanism 
we can not be more precise to say what n* should be. However, if n^ is known 
somehow, the efolding number A^,,, of inflation after horizon crossing can be 
estimated from the assumption that the reheat temperature is of the same 
order as but slightly lower than the energy scale of inflation, Treh ~ 3 x 10^^ 
GeV, through the formula[3]. 



A/", ^ 53 + In 



T 

^ 7 



reh 



1013 GeV 



<57. 



(46) 



We will simply set n* = 2. Then A^^, ^ 54.83847, compatible with the 
constraint from cosmological observation. The pivot wave number of the 
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Fourier mode perturbation is taken as k = a{0)H{0) for which the initial 
efold time nj„j can be defined by a cutoff equation |28j 



' - *<"' -C, (47) 



a{nini)H{nini) e^^^^Hiriini 

The constant Cq should be sufficient large in order that the probable in- 
terferences between the initial curvature and isocurvature modes remain 
negligible [29]. We take Cq ~ 100 for concreteness, which leads to riini = 
—4.60516. To take into account the statistical independence of the initial 
perturbations in our numerical procedure, we integrate Eqs. (l37|) and (138!) 
twice: first with the initial conditions, 



^s{nini) ^ Ss{nini) ^ 



(48) 



and second with the initial conditions. 






(49) 



with Qini = Ssini = l/a{nini)v2k. The curvature and isocurvature pertur- 
bations are customarily described by 7^ = Qcr/<^ and S = 6s/&, for which 
we get TZi and Si [i = 1, 2) after finishing two individual integrations of 
Eqs.f l37|) and fl38l) . T^i (^2) can be thought of as the curvature (isocurvature) 
perturbation without taking into account the coupling to the isocurvature 
(curvature) perturbation, and 7^2 ('^i) does just come from such a coupling. 
The final power spectra and their correlation are calculated as follows [50]: 

^7^ = |^(|7^,p + |7^2^, ^. = |^(|5,p + |52n, 

^7^5 = |^(7^l5l+7^552) . (50) 

In fact, the correlation between the curvature and isocurvature perturbations 
can be simply described by the so-called the relative correlation coefficient 
^ = \%is\/V^Ti^s- The value of "^ lies between and 1, which measures 
to what extent the final curvature perturbation results from the interactions 
with the isocurvature perturbation. 
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Our result of numerical integration for linear perturbations is displayed 
in Figure HI where the power spectra of curvature and isocurvature pertur- 
bations as well as the correlation between them are plotted as functions of 
the number n of efolds. On super-Hubble scales the relative correlation co- 
efficient continues to increase, which is displayed in Figure [51 

The spectra drop off sharply on subhorizon scales until the supposed 
Hubble exit for the Fourier mode perturbation is achieved. On superhori- 
zon scales, however, the spectra become sufficiently flat for 2 < n < 30, 
with S^Ti K. S^g ^ 3.96 X 10~^° and correlator ^7^5 increases steadily from 
1.46 X 10-12 ^o 5 041 X 10-" (3.6 x 10-^ < ^ < 0.132). During the pe- 
riod of 72 > 30, though '^Tis increases at flrst but then drops off steeply, ^ 
increases monotonously. For n > 30, the spectrum of curvature perturba- 
tion begins to increase while the spectrum of isocurvature begins to decrease 
once more. The persistent increase of the curvature perturbation spectrum 
for efold time interval 30 < n < 56.6 reflects a fact that the correlation be- 
tween the curvature and isocurvature perturbations during this time interval 
has become sufficiently strong that it can no longer be neglected. The decay 
of isocurvature perturbation spectrum accelerates in the final stage of infla- 
tion, which becomes unbelievable steep within the last efold and results in 
an unbelievable steep decay of correlator %zs- For n > 56.6, in particular, 
both ^Tis and J^s become negligiblqj, and J^-jz stops varying. Therefore, in 
the present two-fleld model, the power spectrum of curvature perturbation 
on superhorizon scales is changeable with respect to time. This is in con- 
trast with the single fleld inflation case where the curvature power spectrum 
remains approximately a constant after Hubble crossing. Figure [6] gives the 
details of our numerical result for power spectrum of the curvature pertur- 
bation, where the three curves respectively stand for the complete curvature 
spectrum ^-ji = ^ (|7^iP + 17^2^) and its two terms ^|7?.ip and ^|7^2p- 
We see that for 2 < n < 30 on superhorizon scales the spectrum of curvature 
perturbation originates almost exclusively from the fluctuations of the initial 
curvature perturbation along the inflationary trajectory. While for n > 30, 



* If we define a parameter 

a = -^ -=— (51) 

^n + ^s 

to measure the relative contribution of the isocurvature fluctuations, a|n=56. 83847 ~ 3.7 x 
10~^°. Hence the constraint on isocurvature perturbations from observations [31] is trivially 
matched and there is actually no observable isocurvature perturbation spectrum after 
inflation in the present two-field model. 
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the curvature-isocurvature correlation becomes important, the spectrum at 
the final stage of infiation contains considerable contributions from the inter- 
actions with the isocurvature modes. As a result, the curvature perturbation 
spectrum after infiation is slightly higher than the curvature perturbation 
spectrum at Hubble crossing. 

The obtained power spectrum of the curvature perturbation in the present 
two-field model is in good agreement with the COBE normalization of power 
spectrum [32] 

for 2 < n < 30. n ?« 2 is an instant approximately 55 efoldings before the 
end of infiation, it is near the COBE normalization point. So n ~ 2 is a 
reasonable efold time corresponding to Hubble crossing. The scalar spectral 
index of the curvature power spectrum 

n, = l + d\n^Tz/d\nk (53) 

on superhorizon scales is also time-dependent, which has been calculated 
numerically with partial results given in Figure [71 At the claimed Hubble 
crossing, ns\n^2 ~ 0.96, which is compatible with the best observational 
constraint available at present [33] . Ug = 0.98 ± 0.02. 

It follows from the above numerical investigation that, in the two-field 
model under consideration, strong interactions exist between the curvature 
and isocurvature perturbations during the final stage of infiation. Because 
the relative correlation coefficient increases steadily from the assigned tiny 
value deep inside the Hubble radius to values very close to one at the end of 
infiation, a great deal curvature perturbation after infiation originates from 
interactions between the curvature and isocurvature perturbations on super- 
Hubble scales, not only from the quantum fiuctuations of the initial curvature 
perturbation along the classical trajectory. At this point, our model is some- 
what similar to the Double infiation with non-canonical kinetic terms and 
Roulette infiation studied in [23]. What distinguishes our model from the 
Double infiation with non-canonical kinetic terms and Roulette infiation is 
that in the former the power spectrum of isocurvature perturbation remains 
almost a constant during the most of infiationary epoch while in the latter 
it decays rapidly after the Hubble radius crossing. This implies that in our 
model the isocurvature perturbation continues to interact with the curvature 
perturbation on superhorizon scales so that there is no an effective single 
field treatment available for it [3]. 
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4 Discussion 

We have established a multifield infiationary model on the large volume fiux 
compactification scheme in type IIB superstring orientifolds of h^^''^^ > h^^'^\ 
The Kahler moduli are categorized into volume modulus, heavy moduli and 
light moduli. Almost all of closed string moduli emerging from the fiux 
compactification are frozen at the potential minimum in the spirit of KKLT 
mechanism in the large volume limit. However, this is not the case for light 
Kahler moduli. The light Kahler moduli have very small masses, they are not 
fixed by either the 3-form fiuxes or nonperturbative effects in the effective 
superpotential, and could be the driving force for the potential inflationary 
evolution of the universe in its early history. To have a multifleld inflation 
model, the superpotential is expected to contain several exponent terms for 
the light Kahler moduli. Hence we have considered the orientifolds of h^^'^^ > 
4 and supposed that in our constructions the number of light Kahler moduli 
is equal to or greater than two. In some sense, in this paper we have proposed 
a plausible generalization of the better racetrack model P] in the large volume 
approach |71 E]. 

The quantitative analysis of inflationary property in the resulting multi- 
field inflationary model has been done in a simplifled two-fleld model where 
there are only two light Kahler moduli with the axionic components being 
stabilized at the potential minimum also. This assisted-like model is an ex- 
tension of the Kahler moduli inflation model [1] to two-fleld inflation, which 
distinguishes itself by its kinetic Lagrangian that is neither of canonical nor 
non-standard type studied in [23]. Our investigation consists of two steps. 
First, by some inevitable flne tuning in the initial conditions of the scalar 
flelds, we integrate the background equations numerically and obtain the in- 
flationary trajectory of the scalar flelds. As a remarkable characteristic of 
the multifleld model, the trajectory in fleld space in our model are strongly 
bent roughly 50 efolds after the Hubble crossing. The flrst slow-roll param- 
eter is numerically calculated which turns out to almost vanish for the flrst 
60 efolds, implying that the scalar potential is very flat along the trajectory. 
Second, we derive the evolution equations of the linear scalar perturbations. 
These equations are numerically integrated on which we have calculated the 
curvature, isocurvature spectra, and the correlation between them. Our nu- 
merical estimation for the power spectrum of the curvature perturbation and 
the corresponding spectral index is in good agreement with the COBE nor- 
malization and the WMAP observation dataset 1321 531 . In our model, the 
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correlation between the curvature and isocurvature perturbations is gradually 
strong on super-Hubble scales, so that the curvature spectrum approximately 
thirty efolds after the Hubble exit is remarkably different from the curvature 
perturbation spectrum at the Hubble exit, and there might be significant 
non-Gaussianity in bispectrum and trispectrum of the primordial perturba- 
tions. The investigation of the deviations from Gaussianity in power spectra 
in our two-field inflation model is in progress. 
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Figure 2: Classical inflationary trajectory of the two-field model under consid- 
eration, for efolds —4.60516 < n < 56.83847. The inflationary epoch starts 
at n ~ —4.60516 where we put the initial conditions c^mj ~ 5.67 x 10^^ and 
Xini ~ —1.13 X 10~^. In the flrst 35 efolds the trajectory is approximately along 
the ((^ — x)-direction (red curve). When n > 30, the flelds {if it x)/\/2 begin to 
lose their synchrony in evolution. A steep turn in fleld space takes place roughly 
at n PS 55.85 where ip « 5.04 x 10"^, x ~ -2.57 x 10~^ and after that {cp - x)V2 
does no longer evolve but {ip + x)/V^ continues to decrease. During the last efold, 
55.85 < n < 56.83847, the trajectory is along the (93 + x)-direction (magenta 
curve). The inflation approaches to its end at n ~ 56.83847, where if ~ 3.34 x 10~^, 
X ~ -1.94 X 10-3 and ei « 0.97. 
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Figure 3: The first slow-roll parameter ei during the efold time interval 
—4.60516 < n < 56.83847 for our two-field model, ei is much smaller than one for 
most of the efolds, and particularly ei ~ 4.82 x 10~^^ at n = 2, implying that the 
potential is almost flat during the epoch, ei increases drastically during the last 
several efolds. When n = 56.83847, ei ~ 0.97, the inflation tends to end. 



Power Spectra & Correlation of perturbations 
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Figure 4: Power spectra and correlation of perturbations. The red and magenta 
curves represent respectively the spectra of curvature and isocurvature perturba- 
tions. The orange curve describes the correlation between them. All these curves 
are plotted in the efold time interval —4.06516 < n < 56.83847. 
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Correlation Coefficient 
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Figure 5: Relative correlation coefficient '^ in the present model. It vanishes 
initially, as expected. At n = -1, 0, 1 and 2, '(T ps 5.18 x 10"^, 1.35 x 10"^ 1.69 x 
10"'^ and 3.6 x 10~^, respectively. The coefficient increases monotonously on super- 
Hubble scales. When inflation goes close to its end, n k, 56.83847, ^ k, 0.66. 



Power Spectrum of Curvature Perturbation 
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Figure 6: Power spectrum of curvature perturbation in efold time interval 
—4.60516 < n < 56.83847. Red curve describes the evolution of the complete 
curvature spectrum. Pink and magenta curves stand for the partial curvature 
spectra from fluctuations along the inflationary trajectory and from the interac- 
tions with the isocurvature modes, respectively. The curvature spectrum at the end 
of inflation is ^7e|raRi56. 83847 ~ 6.79 x 10"^'', which is different from the curvature 
spectrum =^7^|n~2 ~ 3.96 x lO"^'' at Hubble-exit. 
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Spectral index «, at Hubble crossing 
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Figure 7: Spectral index Ug of the curvature power spectrum in the vicinity of 
Hubble exit. The instant for COBE-like normalization corresponds roughly to 
n = 2, where n^ ~ 0.96. However, the index running is sufficiently small so that 
for 1.8 < n < 50, almost all values of n^ are compatible with the observational 
constraint. 
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